Physica C : Superconductivity XX, XXXXXX (2012) 



Type-1.5 superconductivity in multiband systems: 
magnetic response, broken symmetries and microscopic theory. A brief overview. 

E. Babaevi'2, J. Carlstromi'^, J. Garaud^, M. Silaev^'^ and J.M. Speight"^ 
^Department of Theoretical Physics, The Royal Institute of Technology, Stockholm, SE- 10691 Sweden 
^ Department of Physics, University of Massachusetts Amherst, MA 01003 USA 
^ Institute for Physics of Micro structures RAS, 603950 Nizhny Novgorod, Russia. 
^ School of Mathematics, University of Leeds, Leeds LS2 9JT, UK 



A conventional superconductor is described by a single complex order parameter field which has 
two fundamental length scales, the magnetic field penetration depth A and the coherence length ^. 
Their ratio k determines the response of a superconductor to an external field, sorting them into 
two categories as follows; type-I when n < and type-II when n > 1/V2. We overview here 

multicomponent systems which can possess three or more fundamental length scales and allow a 
separate "type-1.5" superconducting state when, e.g. in two-component case ^1 < V2X < ^2. In 
that state, as a consequence of the extra fundamental length scale, vortices attract one another at 
long range but repel at shorter ranges. As a consequence the system should form an additional 
Semi-Meissner state which properties we discuss below. In that state vortices form clusters in 
low magnetic fields. Inside the cluster one of the component is depleted and the superconductor- 
to-normal interface has negative energy. In contrast the current in second component is mostly 
concentrated on the cluster's boundary, making the energy of this interface positive. Here we 
briefly overview recent developments in Ginzburg-Landau and microscopic descriptions of this state. 
Prepared for the proceedings of Vortex VII conference, Rhodos September 2011. 



Type-I superconductors expel weak magnetic fields, while strong fields give rise to formation of macroscopic normal 
domains with magnetic flux.^^ The response of type-II superconductors is different;"^ below some critical value i^ci, 
the field is expelled. Above this value a superconductor forms a lattice or a liquid of vortices which carry magnetic 
flux through the system. Only at a higher second critical value, Hc2 superconductivity is destroyed. 

These different responses are usually viewed as consequences of the vortex interaction in these systems, the energy 
cost of a boundary between superconducting and normal states and the thermodynamic stability of vortex excita- 
tions. In a type-II superconductor the energy cost of a boundary between the normal and the superconducting state 
is negative, while the interaction between vortices is repulsive.^ This leads to a formation of stable vortex lattices 
and liquids. In type-I superconductors the situation is the opposite; the vortex interaction is attractive (thus mak- 
ing them unstable against collapse into one large "giant" vortex), while the boundary energy between normal and 
superconducting states is positive. 

One can distinguish also a special "zero measure" boundary case where has a critical value exactly at the type- 
I/type-II boundary, which in the most common GL model parameterization corresponds to = 1/a/2- In that case 
vortices do not interact^'^ in the Ginzburg-Landau theory. The noninteracting regime, which is frequently called 
"Bogomolny limit" is a property of Ginzburg-Landau model where, at k = l/\/2, the core-core attractive interaction 
between vortices exactly cancels the current-current repulsive interaction.^'^ However indeed in a realistic system 
even in the limit hz = 1/a/2, there will be always leftover inter- vortex interactions, appearing beyond the GL field 
theoretic description, form underlying microscopic physics. The form of that interaction potential is determined not 
by the fundamental length scales of the GL theory but by non-universal microscopic physics and it can indeed be non- 
monotonic.^ These microscopic corrections are extremely small, however they can be relevant not only at k, = l/\/2 
but also in a very narrow window of parameters near k, ^ 1/ \f2 where intervortex forces in GL theory are also very 
small. We do not consider the physics which arises in the Bogomolny limits in this paper, rather concentrating on 
the physics associated with fundamental modes of GL field theory. 

Recently there has been increased interest in superconductors with several superconducting components. The 
Ginzburg-Landau free energy functional for multicomponent system has the form 



Here are complex superconducting components, I) = V + zeA, and i\)a — |V^a|e*^'', a = 1, 2, and V(|?/^i|) stands for 
effective potential. Depending on symmetry of the system there can also be present mixed (with respect to components 
gradient terms (for a more detailed review see^). 
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FIG. 1. A schematic picture of magnetization curves of type-I, type-II and type-1. 5 superconductors. 



The main situations where multiple superconducting components arise are (i) multiband superconductors'-^^ (where 
tjji represent condensates belonging to different bands), (ii) mixtures of independently conserved condensates such as 
the projected superconductivity in metallic hydrogen and hydrogen rich alloys, ^'^^ (where t/ji represent electronic 
and protonic Cooper pairs or deuteronic condensate) or models of nuclear superconductors in neutron stars interior^^ 
(where i/ji represent protonic and hyperonic condensates) and (iii) superconductors with other than s-wave pairing 
symmetries. The principal difference between the cases (i) and (ii) is the absence of the intercomponent Josephson 
coupling in case of system like metallic hydrogen (ii) because there the condensates are independently conserved. 
Thus the symmetry is U{1) x U{1) or higher. In the case (i) multiple superconducting components originate from 
Cooper pairing in different bands. Because condensates in different bands are not independently conserved there is a 
rather generic presence of intercomponent Josephson coupling |('0iV^2 + '^2i^i) in that case. 



A. Type-1. 5 superconductivity 



The possibility of a new type of superconductivity, distinct from the type-I and type-II in multicomponent 
systems^'^^^' comes from the following considerations. In principle the boundary problem in the Ginzburg-Landau 
type of equations in the presence of phase winding is not reducible to a one-dimensional problem in general. Fur- 
thermore, as discussed in,^'^^'^^'^^^^ in general in two-component models there are three fundamental length scales: 
magnetic field penetration length A and two characteristic length scales of the variations of the density fields , ^2 
which renders the model impossible to parametrize in terms of a single dimensionless parameter Hi and thus the type- 
I/type-II dichotomy is not sufficient for classification. Rather, in a wide range of parameters, as a consequence of the 
existence of three fundamental length scales, there is a separate superconducting regime with ^i/V2 < A < 6/^2- In 
that regime a situation is possible where vortices have long-range attractive (due to "outer cores" overlap) , short-range 
repulsive interaction (driven by current-current and electromagnetic interaction) and form vortex clusters immersed 
in domains of two-component Meissner state. ^^'^^ Recent experimental works^^'^^ proposed that this state is realized 
in the two-band material MgB2. In Ref.^^ this regime was termed "type- 1.5" superconductivity by Moshchalkov et 
al. These works resulted in increasing interest in the subject. ^^^^ Recently type-1. 5 superconductivity was discussed 
in context of quantum Hall effect. 

If the vortices form clusters one cannot use the usual one-dimensional argument concerning the energy of 
superconductor-to-normal state boundary to classify the magnetic response of the system. First of all, the en- 
ergy per vortex in such a case depends on whether a vortex is placed in a cluster or not. Formation of a single isolated 
vortex might be energetically unfavorable, while formation of vortex clusters is favorable, because in a cluster where 
vortices are placed in a minimum of the interaction potential, the energy per flux quantum is smaller than that for 
an isolated vortex. 

Thus, besides the energy of a vortex in a cluster, there appears an additional energy characteristic associated with 
the boundary of a cluster. 

We summarize the basic properties of type-I, type-II and type- 1.5 regimes in the table I.^ 
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single-component Type-I 


single-component Type-II 


multi-component Type- 1.5 


Characteristic 
lengths scales 


Penetration length A & coher- 
ence length e (f < 71) 


Penetration length A & coher- 
ence length e (f > 7f) 


Two characteristic density vari- 
ations length scales ^1,^2 and 
penetration length A, the non- 
monotonic vortex interaction oc- 
curs in these systems in a large 
range of parameters when ^1 < 
V2A<6 


Intervortex inter- 
action 


Attractive 


Repulsive 


Attractive at long range and repul- 
sive at short range 


Energy of su- 
perconduct- 
ing/normal state 
boundary 


Positive 


Negative 


Under quite general conditions 
negative energy of superconduc- 
tor/normal interface inside a vortex 
cluster but positive energy of the 
vortex cluster's boundary 


The magnetic field 
required to form a 
vortex 


Larger than the t her mo dynam- 
ical critical magnetic field 


Smaller than thermodynamical 
critical magnetic field 


In different cases either (i) smaller 
than the thermodynamical critical 
magnetic field or (ii) larger than 
critical magnetic field for single vor- 
tex but smaller than critical mag- 
netic field for a vortex cluster of a 
certain critical size 


Phases in external 
magnetic field 


(i) Meissner state at low fields; 

(ii) Macroscopically large nor- 
mal domains at larger fields. 
First order phase transition be- 
tween superconducting (Meiss- 
ner) and normal states 


(i) Meissner state at low fields, 

(ii) vortex lattices/liquids at 
larger fields. Second or- 
der phase transitions between 
Meissner and vortex states and 
between vortex and normal 
states at the level of mean-field 
theory. 


(i) Meissner state at low fields 

(ii) "Semi-Meissner state": vor- 
tex clusters coexisting with Meiss- 
ner domains at intermediate fields 

(iii) Vortex lattices/liquids at larger 
fields. Vortices form via a first or- 
der phase transition. The transition 
from vortex states to normal state 
is second order. 


rL/uergy Ji/(^i\l ) oi 
N- quantum ax- 
ially symmetric 
vortex solutions 


< jv— 1 Vor- 
tices collapse onto a single N- 
quantum mega- vortex 


N > N-1 ^• 

N-quantum vortex decays into 
N infinitely separated single- 
quantum vortices 


There is a characteristic number Nc 
such that "^^^^ < ^(^^-i) for N < 

Ne, While ^(,^> > for N > 
Nc. N-quantum vortices decay into 
vortex clusters. 



TABLE L Basic characteristics of bulk clean superconductors in type-I, type-II and type-1.5 regimes. Here the most common 
units are used in which the value of the GL parameter which separates type-I and type-II regimes in a single-component theory 
is Kc = Magnetization curves in these regimes are shown on Fig. 1 



B. Generalization to N-component case. 

The concept of type-1.5 superconductivity has a straightforward generalization to N-component case. There it can 
occur in systems where characteristic length scales are ^1, .., ^/e < a/2 A < ^/c+i, ^at and there are thermodynamically 
stable vortices with non-monotonic interaction. 



II. THE TWO-BAND GINZBURG-LANDAU MODEL WITH ARBITRARY INTERBAND 
INTERACTIONS. DEFINITION OF THE LENGTH SCALES AND TYPE-1.5 REGIME 

A. Free energy functional 

In this section we study the type-1.5 regime using the following two-component Ginzburg-Landau (TOOL) free 
energy functional. 
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F = + i(D^2)P^2)* - iyRe{{DxPi){DxP2r} + ^(V x A)^ + Fp (2) 

Here D = V -\- zeA, and ipa = |'^a|e*^'', a = 1,2, represent two superconducting components which, in a two-band 
superconductor are associated with two different bands. The term Fp can contain in our analysis an arbitrary 
cohection of non-gradient terms representing various inter and intra-band interactions. Below we show how three 
characteristic length scales are defined in this two component model (two associated with densities variations and the 
London magnetic field penetration length). Note that existence of two bands in a superconductor is not a sufficient 
conditions for a superconductor to be described by a model like (2) with two well-defined coherence lengths. Conditions 
of appearance of regimes when the system does not allow a description in terms of two-component fields theory (2) is 
discussed in the work based on microscopic considerations. ^^'^^ This kind of two-band GL models were also derived 
earlier from microscopic two-band theories at elevated temperature close, but not too close to Tc.^^^^ Although 
some other recent works proposed to keep higher order gradient terms in GL expansion, however in Ref.^^ it was 
demonstrated that the minimal two-band model (2) is microscopically justified on formal grounds and gives a very 
good description of the system in a wide range of temperatures. 

The only vortex solutions of the model (2) which have finite energy per unit length are the integer A'-flux quantum 
vortices which have the following phase windings along a contour / around the vortex core: V^i = 27rA^, V^2 = 
27rA' which can be denoted as (N,N). Vortices with differing phase windings (N,M) carry a fractional multiple of the 
magnetic flux quantum and have energy divergent with the system size.^^ 

In what follows we investigate only the integer flux vortex solutions which are the energetically cheapest objects 
to produce by means of an external field in a bulk superconductor. Note that since this object is essentially a bound 
state of two vortices, it in general will have two different co-centered cores. 



III. INTERVORTEX FORCES AT LONG RANGE 



In this section we explain how the nature (attractive or repulsive) of the forces between well separated vortices 
in system (2) can be determined purely by analyzing Fp and how three fundamental length scales can be defined in 
the model (2). Below we will analyze system (2) in the case u = but for an arbitrary effective potential. Detailed 
discussion of the effects of mixed gradient terms can be found in^ By gauge invariance, Fp may depend only on 
|?/^2| and S = Oi — We consider the regime that it has a global minimum at some point other than the one with 
IV^al = 0. We may assume, without loss of generality, that the minimum of F is at (iV^il, |'^2|,^) = ('^i,'^2,0) where 
Ui > and U2 > 0. Then the model has a trivial solution, ipi = ui^ ?/^2 = '^2, ^ = 0, which we call the ground state. 
It also supports vortex solutions of the form 

V'a = /a(r)e*^ {A,,A2) = ^{-sme,cos9) (3) 

r 

where /i,/2,« are real profile functions with boundary behavior /a(0) = a(0) = 0, /a(oo) = Ua, a(oo) = — 1/e. No 
explicit expressions for fa, a are known, but, by analyzing the system of differential equations they satisfy, one can 
construct asymptotic expansions for them at large r, see.^'^^ 

At large r from the vortex in the model (2) the system recovers (up to exponentially small corrections) the ground 
state. In fact, the long-range field behavior of a vortex solution can be identified with a solution of the linearization 
of the model about the ground state, in the presence of appropriate point sources at the vortex core. This idea is 
explained in detail for single component GL theory in.^^ A common feature of topological solitons (vortices being a 
particular example) is that the forces they exert on one another coincide asymptotically (at large separation) with 
those between the corresponding point sources interacting via the linearized field theory. For (2), the linearization 
has one vector {A) and 3 real scalar (ei = IV^il — i^i, €2 = \i^2\ — U2 and 5) degrees of freedom. The isolated 
vortex solutions have, by definition within the ansatz we use, S = everywhere. Note that the GL system may also 
possess non-axially-symmetric solutions, such as vortex clusters, and for these there is no reason why S should vanish 
everywhere and in fact it does not.^^ However below we first consider a single vortex using a axially- symmetric ansatz 
, and hence have no source for J, so we can set (5 = in the linearization, which becomes 

Fa„ = ^|Veir + ^|Ve2|^ + ^l • H 1^ ^ \ + hd^A^ - d^A^)' + ^Hul + ul)\A\\ (4) 




Here, H is the Hessian matrix of F^d^/^il, |?/^2|, 0) about ('^1,1^2), that is. 



^ab — 



(5) 



(ni,n2,0) 
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Note that, in Fun, the vector field A decouples from the scalar fields and mediates a repulsive force between vortices 
(originating in current-current and magnetic interaction) with decay length which is the penetration length X = 1/fiA 
where jua is the mass of the field, that is. 



fiA = e^^ul^ul. (6) 

By contrast, the scalar fields ei, 62 are, in general, coupled (i.e. in general the symmetric matrix H has off-diagonal 
terms). To remove these we make a linear redefinition of fields, expanding (61,62)"^ with respect to the orthonormal 
basis for formed by the eigenvectors vi^V2 of 

(^1,62)^ = Xl^l +X2^2. (7) 

The corresponding eigenvalues /ii,/i2 necessarily real (since H is symmetric) and positive (since {ui^U2) is a 
minimum of F^), and hence 
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Fun = IY1 + ^aXa) + ^(^^^^ " S^Al)^ + ^e(u? + ul)\A\^ . (8) 

a=l 

The scalar fields X17X2 each mediate an attractive force between vortices, with length scales 

6 = l/Mi , 6 = 1/M2 (9) 

respectively. Physically these interactions are associated with the attractive core-core interactions. We can be some- 
what more quantitative. In terms of the normal- mode fields X17X2 and A, the composite point source which must be 
introduced into Fun to produce field configurations identical to those of vortex asymptotics is 

Ki=qiS{x), i^2 = q2S{x), J = m{d2, -di)S{x) , (10) 

where hzi is the source for xi, 1^2 the source of X2, j the source for A, d{x) denotes the two dimensional Dirac delta 
distribution and ^1,^2 and m are unknown real constants which can, in principle, be determined numerically by a 
careful analysis of the vortex asymptotics. Physically, a vortex, as seen from a long distance can be thought of as a 
point particle carrying two different types of scalar monopole charge, ^1, <72, inducing fields of mass /ii, /i2 respectively, 
and a magnetic dipole moment m oriented orthogonal to the X1X2 plane, inducing a massive vector field of mass /i^. 
The interaction energy experienced by a pair of point particles carrying these sources, held distance r apart, is easily 
computed in linear field theory. For example, two scalar monopoles of charge q inducing fields of mass ja held at 
positions y and y in experience interaction energy 

Eint = - [ ^x = - [ qS{^ - y)^KMy - yI) = -^^My - y|) (11) 

where n is the source for the monopole at y, x is the scalar field induced by the monopole at y^^ and Kq denotes the 
modified Bessel's function of the second kind. The interaction energy for a pair of magnetic dipoles may be computed 
similarly. In the case of our two component GL model, the total interaction energy has three terms, corresponding to 
the three sources in the composite point source (10), and turns out to be 

222 
Ei^t = "^KoipAT) - ^Koiiiir) - ^Ko{p2r). (12) 

Note that, the first term in this formula which originates in magnetic and current-current interaction is repulsive, 
while the other two as associated with core-core interactions of two kinds of cores are attractive. At very large r, 
Einti^) is dominated by whichever term corresponds to the smallest of the three masses, jua^ A^i, f^2^ so to determine 
whether vortices attract at long range, it is enough to compute just these masses. 

To summarize, the nature of intervortex forces at large separation can be determined purely by analyzing Fpi one 
finds the ground state ('^1,1^2) and the Hessian 1-L of Fp about (1^1,1^2)- From this one computes the mass of the 
vector field A, fiA = ^\/'^ + '^2^ the masses /ii,/i2 of the scalar normal modes (the fields X17X2), these masses 
being the square roots of the eigenvalues of If either (or both) of /ii,/i2 are less than /i^, then the dominant 
interaction at long range is attractive (i.e. vortex core extends beyond the area where magnetic field is localized), 
while if \iA is less than both \i\ and \i2^ the dominant interaction at long range is repulsive. The special feature of 
the two-component model is that the vortices where core extends beyond the magnetic field penetration length are 
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thermodynamically stable in a range of parameters and moreover one can have a repulsive force between the vortices 
at shorter distances. ^'^^'^^ It is important to stress that length scales are not directly associated with the 

individual condensates ?/^2- Rather they are associated with the normal modes defined as^'^^ 

Xi = (IV^il -^^i)cos6- (|V^2| -iX2)sin6, X2 = -(IV^il -^ii)sin6- (|V^2| -iX2)cos6. (13) 

These may be thought of as rotated (in field space) versions of ei = IV^il — i^i, 62 = |'^2| — '^2- The mixing angle^ that 
is, the angle between the x ^^id e axes, is 6, where the eigenvector vi of 1-L is (cos 9, sin 6)^. This, again, can be 
determined directly from 1-L. 

Note also that the shorter of the length scales although being a fundamental length scale of the theory, 

can be masked in a density profile of a vortex solution by nonlinear effects. This, for example certainly happens if 
/i^^ <C /iA = (see short discussion in Ref.^). Also note that in general the minimum of the interaction potential 
will not be located at the London penetration length, because it in general will be also affected by nonlinearities. 



A. Passive band superconductors 

To illustrate the analysis presented above, we consider the simple case of a two band superconductor where one of 
the bands is passive, that is, with a potential of the form 

= -ai I V'l I' + y 1^1 1' + "2 1^/'2 1' - 7(V'i + ^1 V'2) (14) 
where aj,/3i,7 are positive constants. Then Fp is minimized when and ?/^2 have equal phase, and have moduli 



IV^il = ^1 = f 1 + ^) . m =U2 = ^u,. (15) 

y Pi \ aia2 J 0^2 



The mass of the vector field A is 



/iA = e^^i^^^ii^ = eui^l + (16) 
The Hessian matrix of Fp about {ui,U2) is 

It is straightforward to compute explicit expressions for the eigenvalues , /Lt^ of this matrix. These are somewhat 
complicated, but power series expansion in 7 reveals that 

Ml = 27^7+0(7'), M2 = 72^ + 0(72). (18) 

Similarly, the normalized eigenvector associated with eigenvalue is 

^i=( ^ ,-1 ]^0{j') (19) 

y -(2ai - a2) '7 / 

so the normal modes of fiuctuation about the ground state are rotated through a mixing angle 

6 = -{2ai - ^2)-^ + 0(7'). (20) 

In comparison with the uncoupled model (7 = 0) then, we see that, for small coupling 7 the length scales 
l//iyL, l//^i, 1/^2 are unchanged to leading order, but the normal modes with which l//ii,l//i2 are associated 
are mixed to leading order. In particular, there are large regions of parameter space where /i2 < Ma < Mi, so that 
vortices attract at long range, even though the active band, ipi, is naively "type 11" (that is, f3i > e^/4). 
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IV. VORTEX CLUSTERS IN A SEMI-MEISSNER STATE AND NON-PAIRWISE INTERVORTEX 

FORCES. 



A. Model 



In this section, following Ref.^^ we consider in more detail two- component Ginzburg-Landau models, with and 
without Josephson coupling r] which directly couples the two condensates (for treatment of other kinds of interband 
couplings see^). When r] = the condensates are coupled electromagnetically. 



= 1,2 



\{V + ieA)iPif + {2ai+pi\iPif)\i^i\' 



1 



(Vx A)2-r,|Vi||^/'2|cos(^2-^i) 



(21) 



Here again, the gauge covariant derivative is D = \/ -\- ieA, and i/ji = l^/^^le*"^^ are complex fields representing the 
superconducting components. As will be discussed in Sec. VI, calculations based on microscopic two-band Eilenberger 
theory show that the model (21) can be used to study vortex physics in multiband superconductor for a wide range 
or parameters. 

We also discuss importance of complicated non-pairwise forces between superconducting vortices arising in certain 
cases in multicomponent systems. ^^'^^ These non-pairwise forces in certain cases have important consequences for 
vortex clusters formation in the type- 1.5 regime. 



B. Vortex clusters in a Semi-Meissner state and non-pairwise interactions. 



In this section we allow fiuctuations in phase difference. When t here is non-zero i nterband Josephson coupli 
phase difference is associated with a massive mode. Its mass is ^/r]{u1 + U2)/uiU2. 



iing. The 




FIG. 2. Ground state of Ny = 9 flux quanta in a type- 1.5 superconductor enjoying U{1) x U{1) symmetry of the potential 
{i.e. Tj = 0). The parameters of the potential being here (ai,/3i) = ( — 1.00, 1.00) and (0^2, A) = (—0.60, 1.00), while the electric 
charge is e = 1.48. The displayed physical quantities are a the magnetic flux density, b (resp. c) is the density of the first 
{resp. second) condensate |'0i,2|^- d (resp. e) shows the norm of the supercurrent in the first (resp. second) component. Panel 
f is Im('0i'02) = |'0i||'02| sin(^2 — ^1) being nonzero when there appears a difference between the two condensates. Parameters 
are chosen so that the second component has a type-I like behavior while the first one tends to from well separated vortices. 
The density of the second band is depleted in the vortex cluster and its current is mostly concentrated on the boundary of the 
cluster (see Ref.^°). 
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Fig. 2 and Fig. 3 show numerical solutions for N-vortex bound states in several regimes (for technical details see 
Appendix of^^). A Nonlinear Conjugate Gradient scheme was used to solve the finite element formulation of the 
variational problem, within the framework provided by the Freefem++ library. Animations showing the evolution 
of the system, during the numerical energy minimization, from the various initial configurations to the vortex clusters 
in the energy minimization process can be found online. The common aspect of the regimes shown on these figures 
is that the density of one of the components is depleted in the vortex cluster and has its current most concentrated 
on the boundary of the vortex cluster (i.e. has a type-I like behavior). At the same time the second component forms 
a full-fiedged vortex lattice inside the vortex cluster (i.e. has a type-II-like behavior). 




FIG. 3. Elongated ground state cluster of 18 vortices in a superconductor with two active bands. Parameters of the interacting 
potential are (ai,/3i) = ( — 1.00, 1.00), (a2,/32) = (—0.0625,0.25) while the interband coupling is 77 = 0.5. The electric charge, 
parameterizing the penetration depth of the magnetic field, is e = 1.30 so that the well in the nonmonotonic interacting 
potential is very small. In this case there is visible admixture of the current of second component in vortices inside the cluster, 
though its current is predominantly concentrated on the boundary of the cluster. 



C. Non-compact vortex clusters and non-pairwise interactions 

Next we report the regime where the passive second band (i.e. with positive 0^2) is coupled to the first band by 
extremely strong Josephson coupling r] = 7.0 (shown on Fig. 4). This coupling imposes a strong energy penalty both 
for disparities of the condensates variations and for the difference between phases of the condensates. Besides that in 
that regime there is a relatively strong non-pairwise interactions which diminish the energetic benefits of a triangle-like 
states compared to line-like vortex states. We get a fiat and complicated energy landscape and the outcome of the 
energy minimization strongly depends on initial configuration. Simulations whose outcome is compact clusters like 
Fig. 2 and Fig. 3 clearly ground states, since various initial guesses lead to similar final configurations. Simulating 
systems like in Fig. 4 is less straightforward. Numerical evolution in these systems is extremely slow because of the 
very complicated energy landscape. The final state strongly depends on the initial field configuration, indicating 
the configuration is not ground state but a bound state with a very slow evolution. Formation of highly disordered 
states and vortex chains due to short-range nature of the attractive potentials and many-body forces was a common 
outcome of the simulation in the similar type- 1.5 regimes with strong Josephson coupling, in spite of negligible effects 
of ultra-fine numerical grid. 

The Fig. 4 shows the typical non-universal outcome of the energy minimization in this case. Striking feature here 
is formation of vortex stripe-like configuration. Indeed it strongly contradicts the ground state expected from the 
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two-body forces in this system. Namely the axiahy symmetric two-body potentials with long range attraction and 
short-range repulsion (which we have in this case) do not allow stripe formation in the ground state configurations. 
Nonetheless such structures are physically entirely possible: first the exponentially decaying short-range forces result 
in very tiny forces which can drive evolution from a stripe to compact cluster configuration. Secondly there are 
repulsive multibody intervortex forces which compromise system's ability for form a compact configuration.^^ So in 
that regime vortex stripes and small lines can easily form because of non-pairwise intervortex interactions. 

Note that even in this regime, the system exhibits self-induced gradients of the phase difference, in spite of the 
strong Josephson coupling. 
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FIG. 4. A bound state of an Nv = 25 vortex configuration in case when superconductivity in the second band is due to 
interband proximity effect and the Josephson coupling is very strong 77 = 7.0. The initial configuration in this simulation was a 
giant vortex. Other parameters are = ( — 1.00, 1.00), (a2, A) = (3.00,0.50), e = 1.30. For the simulations, like the one 

shown on Fig. 3 the stopping criterion of energy minimization was when relative variation of the norm of the gradient of the GL 
functional with respect to all degrees of freedom to be less than 10~^. Here the situation is slightly different from that shown 
on two previous figures. Clearly in the shown above configuration the ground state was not reached. However the interaction 
potentials a such that the evolution at the later stages becomes extremely slow. The number of energy minimization steps in 
this case was order of magnitude larger than what was required for convergence in the previous regimes. 



V. SEMI-MEISSNER STATE AND BROKEN SYMMETRIES. 



As discussed above a system with non-monotonic intervortex interaction potentials allow a state with macroscopic 
phase separation in vortex droplets and Meissner domains. In type-1.5 superconductors this state can also represent 
a phase separation into domains of states with different broken symmetries. In this section we will give two different 
examples of how such behavior can arise. 

Note that in multicomponent superconductors some symmetries are global (i.e. associated with the degrees of 
freedom decoupled from vector potential) and some are local i.e. associated with the degrees of freedom coupled to 
vector potential. It is well known that in the later case the concept of spontaneous symmetry breakdown is not defined 
the same way as in a system with global symmetry. However below, for brevity we will not be making terminological 
distinctions between local and global symmetries (detailed discussion of these aspects can be found in e.g.^^'^^'^^). 
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A. Semi-Meissner state as a macroscopic phase separation in U{1) x U(l) and U(l) domains. 

Consider a superconductor with broken U{1) x U{1) symmetry, i.e. a cohection of independently conserved conden- 
sates with no intercomponent Josephson couphng. As discussed above, in the Semi-Meissner state in vortex droplets 
the superconducting component which has vortices with larger cores is more depleted. In U{1) x U{1) system the 
vortices with phase windings in different condensates are bound electromagnetically which gives asymptotically loga- 
rithmic interaction potential with prefactor proportional to |'0ip|'02p/(|'0iP + |'02p),^^ and even weaker interaction 
strength at shorter separations. 

Consider now a macroscopically large vortex domain. Even if the second component there is not completely 
depleted, its density is suppressed and as a consequence the binding energy between vortices with different phase 
windings (A^i = 27r, AO2 = 0) and (A^i = 0, AO2 = 27r) can be arbitrarily small. Moreover the vortex ordering 
energy in the component with more depleted density will also small. As a result, even small thermal fluctuation can 
drive vortex sublattice melting transition^^'^^in the macroscopically large vortex droplet. In that case the fractional 
vortices in weaker component tear themselves off the fractional vortices in strong component and form a disordered 
state. Note that the vortex sublattice melting is associated with the phase transition from U{1) x U{1) to U{1) 
state. ^^'^^ I.e. that vortex cluster will represent a domain of U{1) phase (associated with the superconducting state 
of strong component) immersed in domain of vortexless U{1) x U{1) Meissner state. If the magnetic field is increased 
the system will go from the Semi-Meissner state (with coexisting U{1) x U{1) and U{1) domains) to U{1) vortex state. 

B. Semi-Meissner state as a macroscopic phase separation in U{1) and U{1) x Z2 domains in three band 

type-1.5 superconductors. 

In the previous subsection we considered the case where Semi-Meissner state represents coexistences of domains 
with different broken symmetries as a consequence of vortex sublattice melting transition. In type-1.5 systems the 
coexistence of domains with different broken symmetry can also take place also in the ground state, i.e. without 
the need of thermal fluctuations.^^ In this subsection we discuss an example studied in Ref.^^ of such behavior in a 
three-band model with "phase frustration" . 

The minimal GL free energy functional to model a three-band superconductor is 

F=^(VxA)2+ J2 ^|DV'i|'+«i|V'i|' + ^/?i|V'i|'+ E T.'^iM'^jloosi^ij). (22) 

i=l,2,3 i=l,2,3 j>i 

Here the phase difference between two condensates are denoted cpij = cpj — cpi. 

Systems with more than two Josephson-coupled bands can exhibit phase frustration. ^^^'^'^''^^ For r]ij < 0, a given 
Josephson interaction energy term is minimal for zero phase difference (we then refer to the coupling as "phase- 
locking" ), while when rjij > it is minimal for a phase difference equal to tt (we then refer to the coupling as 
"phase-antilocking" ). Two component systems are symmetric with respect to the sign change rjij —r]ij as the 
phase difference changes by a factor tt, for the system to recover the same interaction. However, in systems with more 
than two bands there is generally no such symmetry. For example if a three band system has 77 > for all Josephson 
interactions, then these terms can not be simultaneously minimized, as this would correspond to all possible phase 
differences being equal to tt. 

The ground state values of the fields and (pij of system (22) are found by minimizing its potential energy 

} -^^Vij\i^i\\i^j\cos{(pij). (23) 

i j>i 

Minimizing the potential energy (23) can not in general be done analytically. Yet, some properties can be derived 
from qualitative arguments. In terms of the sign of the ?7's, there are four principal situations: 



Case 


Sign of 7^12, ^13, ^23 


Ground State Phases 


1 




^1= ^2 = ^3 


2 


- - + 


Frustrated 


3 


- + + 


V^l = V^2 = ^3 + TT 


4 


+ + + 


Frustrated 



The case 2) can result in several ground states. If \r]23\ <C |^i2|, 1^131, then the phase differences are generally 
(Pij = 0. If on the other hand \r]i2\, \Vi3\ 1^23 1 then (p23 = tt and (fi2 is either or tt. For certain parameter values 
it can also have compromise states with cpij not being integer multiples of tt. 
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The case 4) can give a wide range of ground states, as can be seen in Fig. 5. As rji2 is scaled, ground state phases 
change from (— tt, tt, 0) to the hmit where one band is depleted and the remaining phases are (— 7r/2, 7r/2). 



Ground state phases 













— w 




a) b) 
- - -O O- ■ 




2.5 



3.5 



4.5 



'12 



FIG. 5. Ground state phases of the three components as function of 7712 (here (/?3 = fixes the gauge). The GL parameters 
are = 1, /^^ = 1, 7713 = 7723 = 3. For intermediate values of 7712 the ground state exhibits discrete degeneracy (symmetry 
is U{1) X Z2 rather than U{1)) since the energy is invariant under the sign change (p2 ^ —^2, ^3 ^ — ^3. For large 7712 we 
get (p2 — ^3 = TV implying that {ipsl =0 and so there is a second transition from U{1) x Z2 to U{1) and only two bands at the 
point d). Here, the phases were computed in a system with only passive bands, though systems with active bands exhibit the 
same qualitative properties except for the transition to U{1) and two bands only {i.e. active bands have non-zero density in 
the ground state). 

An important property of the potential energy (23) is that if any of the phase differences cpij is not an integer multiple 
of TT, then the ground state posses an additional discrete Z2 degeneracy. For example for a system with = — 1, ft = 1 
and r]ij = 1, two possible ground state are given by (pi2 = 27r/3, (pis = — 27r/3 or (pi2 = — 27r/3, (pis = 27r/3. Thus in 
this case, the symmetry is U{1) x Z2, as opposed to U{1). As a result, like any other system with Z2 degeneracy, the 
theory allows an additional set of topological excitations : domain walls interpolating between the two inequivalent 
ground states. 

The ground state of a phase frustrated superconductor is in many cases non-trivial, with phase differences being 
compromises between the various interaction terms. Inserting vortices in such a system can shift the balance between 
different competing couplings, since vortices can in general have different effects on the different bands. In particular, 
since the core sizes of vortices are not generally the same in all bands, vortex matter will typically deplete some 
components more than others and thus can alter the preferred values of the phase difference. So the minimal potential 
energy inside a vortex lattice or cluster may correspond to a different set of phase differences than in the vortexless 
ground state. To see this, consider the following argument: The phase-dependent potential terms in the free energy 
(22) are of the form 



(24) 



where Ui are ground state densities and each fi{r) represent an Ansatz which models how superfluid densities are 
modulated due to vortices. Consider now a system where N vortices are uniformly distributed in a domain ft. The 
phase dependent part of the free energy is 



i>j 



/ dTfi{T)fj{T)cOs{(pij{T)). 

Jn 



(25) 



If (pij is varying slowly in comparison with the inter vortex distance, then it can be considered constant in a uniform 
distribution of vortices (as a first approximation). In that case (25) can be approximated by 



(26) 
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If on the other hand Lpij varies rapidly, then it is not possible to define fjij without a spatial dependence. Then ^pij 
will depend on %(r) which is related to the local modulation functions fifj and vary with a characteristic length 
scale. 

Thus, fj is the effective inter-band interaction coupling resulting from density modulation. Since in general, fi ^ fj 
(unless the two bands i, j are identical), one must take into account the modulation functions fi when calculating 
the phase differences. In particular, if the core size in component i is larger than in component j, then J drfifk < 
J drfjfk and therefore the phase differences cpij minimizing (26) depend on and consequently on the density of 
vortices. Roughly speaking, introducing vortices in the system is equivalent to relative effective decrease of some of 
the Josephson coupling constants. 

This can have profound consequences, as the symmetry of the problem depends on the Josephson interaction terms. 
In Figs. 6 , 7 we see a type-1.5 system in which the symmetry of the ground state is U{1). As vortices are inserted 
into the system, they form clusters and the effective inter band interactions fjij are renormalized to a degree that the 
symmetry of the domain near vortex clusters changes to U{1) x Z2. Thus the Semi-Meissner state in such a system 
represents macroscopic phase separation in domains of broken U{1) and U{1) x Z2 symmetries. 



VI. MICROSCOPIC THEORY OF TYPE-1.5 SUPERCONDUCTIVITY 

The phenomenological Ginzburg-Landau model described above predict the possibility of 1.5 superconducting state. 
Formally the GL theory applies only at elevated temperatures. To describe type-1.5 superconductivity in all temper- 
ature regimes (except, indeed the region where mean-field theory is inapplicable) as well as to make a connection with 
a certain class of the real systems requires a microscopic approach which also does not rely on a GL expansion. Such 
a theory was recently developed in.^^ So let us consider the described above physics in a microscopic formalism of 
self-consistent Eilenberger theory. We consider a superconductor with two overlapping bands at the Fermi level. ^ The 
corresponding two sheets of the Fermi surface are assumed to be cylindrical. Within quasi-classical approximation 
the band parameters characterizing the two different sheets of the Fermi surface are the Fermi velocities Vpj and the 
partial densities of states (DOS) vj^ labeled by the band index j = 1,2. We normalize the energies to the critical 
temperature Tc and length to tq = hVpi/Tc. The system of Eilenberger equations for two bands is 

VFjU^ (V + iA) fj + 2ujnfj - 2Ajgj = 0, (27) 
VFjU^ (V - iA) /+ - 2ujnf^ + 2A*^, = 0. 

Here ujn = (2n + l)7rT are Matsubara frequencies and vfj = Vfj/VfI' The vector np = (cosOp^sinOp) parameterizes 
the position on 2D cylindrical Fermi surfaces. The quasi-classical Green's functions in each band obey normalization 
condition + /j// = 1- The Eilenberger differential Eqs.(27) are solved together with the integral self-consistency 
equations for the gaps 

A, = tV / X^jfjdOp (28) 
^ = ^j^Fj'^Im / n^gjdOp. (29) 

^ i = l,2 n=0 

Here Xij is the coupling matrix which satisfies the symmetry relations niAi2 = 1^2X21 where rii are the partial DOS 
normalized so that ni + 77,2 = 1. We consider An > A22 and therefore refer to the first band as "strong" and to the 
second as "weak". 

The asymptotics of the gap functions |Ai^2|(^) at distances far from the vortex core can be found linearizing 
the Eilenberger Eqs.(27) together with the self-consistency equations. The asymptotics of the linearized system is 
governed by the singularities of response function found in Ref.(^^) among which are the poles and branch cuts. In 
general there are two regimes regulating the asymptotic behavior of gap functions. The first regime is realized when 
two poles of the response function lie below the branch cut. The two poles determine the two inverse length scales or, 
equivalently, the two masses of composite gap functions fields (i.e. linear combinations of the fields as in the previous 
section), which we denote as "heavy" jhh and "light" jul (i.e. jj^h > I-^l)- At elevated temperatures these masses 
are exactly the same as the given by a corresponding two-component GL theory obtained by the gradient expansion 
from the microscopic theory. The GL theory can be also used to describe the two massive modes at relatively low 
temperatures.^^ In this case the coefficients should be adjusted phenomenologically in order to obtain the same values 
of masses as the microscopic theory yields. 
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FIG. 6. Interacting vortex clusters with internal Z2 symmetry in a frustrated three band superconductor. The snapshot 
represents a non-stationary state of the weakly interacting well- separated clusters. In this numerical computation, each of 
the clusters has with a good accuracy converged to a physical solution of GL equations, but the snapshot is taken during 
the slow evolution driven by the weak intercluster interaction. The snapshot demonstrates the existence of long-range field 
variations associated with the soft mode. This produces long-range weak intervortex forces. Displayed quantities are: a) 
Magnetic field, b-d) iV'i |'02|^, iV'sl^, e) | |V^2| sin (^12, f) iV^iUV^sl sin(^i3), g) iV^iUV^sl sin(^23). The GL parameters are 
a\ — —3, /3i =3, a2 — —3, fi^ = 3, as = 2, /^s — 0.5, 7712 — 2.25, 7713 — —3.7. The parameter set was chosen so that 
it lies in the regime where the ground state symmetry of the system without vortices is but is close to the x Zi 

region. Because of the disparity in vortex core size the effective interaction strengths f]ij are depleted to different extents. As 
a consequence, the symmetry associated with the effective couplings inside the cluster correspond to the x Z2 regime. 



The second regime is realized at lower temperatures when there is only one pole below the branch cut, In this case 
the asymptotic is determined by the light mass mode /i^ and the contribution of the branch cut which has all the 
length scales smaller than some threshold one determined by the position of the lowest branch cut on the imaginary 
axis. The branch cut contribution is essentially non-local effect which is not captured by GL theory therefore one can 
expect growing discrepancies between effective GL solution and the result of microscopic theory at low temperatures. 

The examples from ref.^^ of the temperature dependencies of the masses /iL,//(^) are shown in the Fig. 8. The 
evolution of the masses \il,h is shown in the sequence of plots Fig.8(a)-(d) for Aj increasing from the small values 
A J <C Aii,A22 to the values comparable to intraband coupling Aj ~ Aii,A22. The two massive modes coexist at the 
temperature interval < T < Tc, where the temperature is determined by the branch cut position, shown in 
the Fig. 8 by black dashed line. For temperatures T < there exists only one massive mode lying below the branch 
cut. At very low temperatures the mass /iL is very close to the branch cut. As the interband coupling parameter is 
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FIG. 7. Interacting vortex clusters with internal Zi symmetry in a frustrated three band superconductor. The panel A 
displays the magnetic field B. Panels B and C respectively display sin(/?i2 and sin(/?i3, the third phase difference can obviously 
be obtained from these two ones. Second line, shows the densities of the different condensates |'0i|^ (D), |'02|^ (E), iV^sl^ (F). 
The third line displays the supercurrent densities associated with each condensate |Ji| (G), IJ2I (H), IJ3I (I). The parameter 
set here is the same as in Fig. 6. The instructive difference here is that th sine of the phase differences is represented 'unweighted' 
by the densities in contrast to Fig. 6. Panel C now makes clear that the inner cluster is in a defined state (/?i3 ~ 7r/2 (whose 
opposite state would have been — 7r/2). Panel B gives a visualization of the long range interaction between the clusters. 



increased, the temperature T^* rises and becomes equal to at some critical value of Aj = Ajc- 

Besides justifying the predictions of phenomenological two-component GL theory^^ the microscopic formalism 
developed in Ref.(^^) allows to describe type-1.5 superconductivity beyond the validity of GL models. The type-1.5 
behavior requires a density mode with low mass /i^ to mediate intervortex attraction at large separations, which 
should coexist with short-range repulsion. In ref.^^ find that the temperature dependence of is characterized 

by an anomalous behavior, which is in strong contrast to temperature dependence of the mass of the gap mode in 
single-band theories. As shown on Fig. 9a the function fiL{T) is non-monotonic at low temperatures. 

Therefore for a certain range of parameters in contrast with the physics of singe-band superconductors the product 
of London penetration depth A and /i^ has a strong and nonmonotonic temperature dependence shown in Fig. 9b. The 
inverse of the mass of the light composite gap mode jj^l sets the range of the attractive density-density contribution 
to intervortex interaction. Therefore the condition for the occurrence of the intervortex attraction will be met if 
A/Jjl < 1. Note that only infinitesimally close to Tc, the product can be interpreted as single-component-like GL 
parameter k, because the inverse mass V^j^J^^ becomes the GL coherence length. However as shown on Fig. 9 b) this 
GL parameter can have a very strong temperature dependence, thus in general it cannot be used as a single parameter 
to characterize two-band systems. This is because in a two-band superconductor, for a wide range of parameters even 
slightly away from Tc the temperature dependence of /i^, is dramatically different from that of the inverse magnetic 
field penetration length A~^. 

Furthermore because the softest mode with the mass /iL in two band system may be associated with only a fraction 
of the total condensate (as follows from corresponding mixing angles), and because there could be the second mixed 
gap mode with larger mass /j^h^ the short-range intervortex interaction can be repulsive. In Ref.(^^) the temperature 
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FIG. 8. Calculated in^^ masses fiL and fiH (red solid lines) of the composite gap function fields for the different values of 
interband Josephson coupling Aj and 7f = 1- In the sequence of plots (a)-(d) the transformation of masses is shown for Aj 
decreasing from the small values Aj ^ An, A22 to the values comparable to intraband coupling Aj ^ An, A22- The particular 
values of coupling constants are An = 0.25, A22 = 0.213 and Aj = 0.0005; 0.0025; 0.025; A22 for plots (a-d) correspondingly. 
By black dash-dotted lines the branch cuts are shown. In (a) with blue dash-dotted lines the masses of modes are shown for 
the case of Aj = 0. Note that at Aj = the two masses go to zero at two different temperatures. Because 1//jl,h are related to 
the coherence length, this reflects the fact that for U{1) x U{1) theory there are two independently diverging coherence lengths. 
Note that for finite values of interband coupling only one mass fiL goes to zero at one Tci this is in turn a consequence of the 
fact that Josephson coupling breaks the symmetry down to single U{1). 



dependencies of and /i^ are compared demonstrating how in these cases the system goes from type-II to type- 1.5 
behavior as temperature is decreased and fiL becomes smaller than A~^, and, the density associated with the light 
mode is small enough that the system has a short-range intervortex repulsion. 

We calculate self-consistently the structure of isolated vortex for different values of = vf2/vfi. A complex 
aspect of the vortex structure in two-band system is that in general the exponential law of the asymptotic behavior 
of the gaps is not directly related to the "core size" at which gaps recover most of their ground state values. We can 
characterize this effect by defining a "healing" length L/\i of the gap function as follows \Ai\{L/\i) = 0.95A^o. The 
characteristic example of the vortex structure is shown in Fig. 9(c). For this case we obtain that Lai ^ 0.8 for all 
values of jf- On the contrary, the healing length I/A2 of changes significantly such that I/A2 = 1-6; 2.5; 3.2; 3.9; 4.5 
for = 1; 2; 3; 4; 5 correspondingly. 

To demonstrate the type-1.5 superconductivity i.e. large-scale attraction and small-scale repulsion of vortices which 
originates from disparity of the variations of two gaps, the intervortex interaction energy was calculated in.^^ The 
two-band generalization of the Eilenberger expression was evaluated for the free energy of the two vortices positioned 
at the points tr = ((i/2,0) and tl = (— (i/2,0) in xy plane. In Fig. 9(d) the interaction energy Eint is shown as a 
function of the distance between two vortices d. The energy Eint is normalized to the single vortex energy Ey. The 
plots on Fig. 9(d) clearly demonstrate the emergence of type-1.5 behavior when the parameter 7i?, which characterizes 
the disparity in band characteristics is increased. This is manifested in the appearance non-monotonic behavior of 
the intervortex interaction energy Eint{d) as a consequence of two-component structure of the theory. 
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FIG. 9. Calculated in ref.^^ (a) Masses fiL and /j^h (red solid and dotted lines) of the composite gap function fields and inverse 
London penetration (blue dashed lines) for the different values of A/j.l{Tc)/\/2 = 1; 2; 3; 5. The position of branch cut is shown 
by black dash-dotted line, (b) The temperature dependence of the quantity A/jl for A/j.l{Tc)/ = 1;2;3;5 (red solid, blue 
dashed and black dash-dotted lines), (c) Distributions of magnetic field H{r)/H{r — 0), gap functions |Ai|(r)/Aio (dashed 
lines) and |A2|(r)/A2o (solid lines) for the coupling parameters An = 0.25, A22 = 0.213 and A21 = 0.0025 and different values 
of the band parameter — 1;2;3;4;5. (d) The energy of interaction between two vortices normalized to the single vortex 
energy as function of the intervortex distance d. In panels (c,d) the temperature is T = 0.6. 



VII. CONCLUSION 

We reviewed the recent developments in description of type- 1.5 superconductivity in multicomponent systems. Both 
at the levels of microscopic and Ginzburg-Landau field theories the behavior arises as a consequence of one or several 
of the fundamental length scales associated with density variations being larger than the magnetic field penetration 
length A, while at the same time the system possessing thermodynamically stable vortex excitations. These vortices 
have long-range attractive (originating from outer cores overlap) and short-range repulsive interaction. This leads to 
an additional "semi-Meissner" phase sandwiched between Meissner and vortex states which is a macroscopic phase 
separation into domains of Meissner and vortex states. We discussed that in case of thermal fluctuations or/and more 
that two bands this phase separation can also results in coexistence of macroscopically large domains with different 
broken symmetries: i.e. U{1) and U{1) x U{1) or U{1) and U{1) x Z2. 
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